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We have explored Pauli paramagnetism. Landau diamagnetism and de Haas-van Alphen effect in 
a single framework, and unified these three effects for all temperatures as well as for all strengths of 
magnetic field. Our result goes beyond Pauli-Landau result on the magnetism of the 3-D ideal gas of 
electrons, and is able to describe crossover of the de Haas- van Alphen oscillation to the saturation of 
magnetization. We also have obtained a novel asymptotic series expansion for the low temperature 
properties of the system. 
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I. INTRODUCTION 

Pauli paramagnetism Landau diamagnetism Q 
and de Haas- van Alphen effect Q are very common topics 
of condensed matter and statistical physics not only for 
the undergraduate and graduate students 0-0| but also 
for the theoreticians [8l-[l3| and experimentalists fl5l - l20| . 
The simplest system of interest, for these particular top- 
ics, is a 3-D ideal gas of electrons exposed in a constant 
magnetic field (B = Bk) 0, 0]. For the existence of 
spin, each electron behaves as a tiny magnet; and being 
the value of spin to be 1/2, the gas of electrons obeys 
Fermi-Dirac statistics. Response of the external mag- 
netic field (Bk) to the tiny magnets (electrons) of dipole 
moment fiss is known as Pauli paramagnetism. On the 
other hand, magnetic field induces orbital motions to the 
charged electrons, and closed electric circuits in the x — y 
plane are produced every where in space obeying Fermi- 
Dirac statistics. These circuits according to the Lenz's 
law oppose the external magnetic flux to pass through 
them, and result Landau diamagnetism. For the case of 
strong magnetic field, apart from Pauli paramagnetism 
and Landau diamagnetism; oscillations of magnetization 
of the free electron gas with period ^ 1/B are observed. 
Such a phenomenon is called de Haas- van Alphen effect. 
These are the inclusive scenario for the magnetization of 
a 3-D ideal gas of electrons, and are very important in 
statistical and condensed matter physics, in particular, 
for the characterization of metals. 

That electrons have nonzero magnetic moment for 
their intrinsic spin (1/2), and that they obey Fermi-Dirac 
statistics, was experimentally verified after Pauli's pre- 
diction of the paramagnetism of the electron gas in a 
metal exposed in a weak magnetic field [l| . Discrepancies 
of Pauli's prediction and experimental data 2lJ for weak 
magnetic field was attributed by Landau's prediction on 
the diamagnetic contribution from the orbital motion of 
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FIG. 1: Thick and thin solid lines in (a) represent Pauli- 
Landau (P-L) results on magnetization of the ideal gas of 
electrons for T /Tf ~ 0.01 and 0.1 respectively. Dashed line 
represents weak field result for T 0. Thick and thin solid 
lines in (b) represent P-L results for ^bB /kBTp = 2.5 and 
0.05 respectively. 



electrons Landau showed, that, this diamagnetic con- 
tribution, to the lowest order in magnetic field, is 1/3 of 
the paramagnetic contribution predicted by Pauli. Thus, 
in presence of a weak magnetic field, ideal gas of electron 
is a paramagnet. Although it was a natural trend to con- 
sider the magnetization of the electron gas as a separate 
sum of Pauli's and Landau's results Landau-Lifshitz 
unified Pauli paramagnetism and Landau diamagnetism 
in a single framework to justify the separate addition of 
the two contributions [6]. But, this unification was made 
only to the lowest order in magnetic field. For strong 
magnetic field. Landau obtained a formula only for the 
oscillatory part of the magnetization, and predicted it to 
be the theory for the de Haas- van Alphen effect. But, he 
did not unify the three effects altogether [255. I* ^ 
natural question, whether the 1/3 relationship valid for 
all strengths of the magnetic field; if not, how it changes 
for all values of the magnetic field. This question was 
partially answered by Sondheimer- Wilson [2^ . They 
unified the three effects in the low temperature regime 
for weak as well as for strong magnetic field, and ob- 
tained the magnetization as a sum of Pauli's weak field 
formula for the paramagnetism. Landau's weak field for- 
mula for the diamagnetism, and Landau's strong field 
formula for the de Haas- van Alphen effect. Although no 
new result, what Landau obtained after the calculation 
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of the de Haas- van Alphen effect [25j, was obtained by 
Sondheimer- Wilson, yet it is a unification of the three ef- 
fects for weak and strong magnetic fields, and it increases 
the range of applicability of Pauli-Landau results for the 
three effects. For this reason, here-from, we represent 
Pauli-Landau result on the 3-D ideal gas of electrons by 
Sondheimer- Wilson formula 23]. 

Being the topic old, before going into the details, we 
should explain our motivation of writing this paper. Let 
us start from Pauli-Landau result (or Sondheimer- Wilson 
formula) as presented in FIG. 1. Sondheimer- Wilson for- 
mula represents partial unification of the three effects ac- 
cording to their approximation scheme. Incompleteness 
of Pauli-Landau result (or Sondheimer- Wilson formula) 
is clearly evident in this figure. It neither explains the 
saturation of magnetization of the ideal gas of electron 
nor it does not have a classical limit. So, Sondheimer- 
Wilson formula needs to be completed for all tempera- 
tures as well as for all strengths of the magnetic field. 
In the following, we will explore Pauli paramagnetism. 
Landau diamagnetism and de Haas- van Alphen effect of 
the 3-D ideal gas of electrons in a single framework, and 
unify these three effects for all temperature we well as for 
all strengths of magnetic field. 

Calculations of this paper will begin with a general- 
ized form of the Landau levels for the 3-D ideal gas of 
electrons exposed in a constant magnetic field. Then 
we will proceed adopting the steps of Landau-Lifshitz 
for unifying Pauli paramagnetism and Landau diamag- 
netism and that for the quantitative introduction of de 
Haas- van Alphen effect Q. Then we will proceed in our 
own way to unify these three effects, and to get the most 
general form of the grand potential of the system. There- 
from we will evaluate magnetization, susceptibility and 
thermodynamic energy. We will plot field and tempera- 
ture dependences of these thermodynamic variables, and 
compare our results with that of Pauli-Landau [23| . 
For the low temperature regime, we will also generalize 
Pauli-Landau result or Sondheimer- Wilson formula in a 
novel way @, [2^ . 



II. GRAND POTENTIAL 
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A. Landau-Lifshitz integral form of grand potential 

Grand potential of our system, considering Landau 
diamagnetic and Pauli paramagnetic contributions, is 
thus given by 
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where /(/i) is given by 
ksTrnV 
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/i is the chemical potential of the system in absence of the 
magnetic field, V = AL is volume of the system, A is area 
along x — y plane of the system, and L is the extent of the 
system along the z direction. Summation in Eqn. ^ can 
be exactly obtained using Poisson summation formula: 
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both sides of this formula by g{x) and integrating over x 
from to oo, we can have 
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Substituting g{j) as /(/i — 2(1 bBj), and putting right 
hand side of Eqn.® into the square bracket in Eqn.Q 
we get 
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Let us consider a system of 3-D ideal gas of electrons 
exposed in a constant magnetic field (B = Bk) to be in 
equilibrium with a heat and particle reservoir of temper- 
ature T and chemical potential /x. Let the average total 
number of particles (electrons) be N, and the mass of 
each particle be m. Energy levels of each particle of the 
system are given by the generalized form of the Landau 
levels |6| 
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where is momentum of the particle along the z direc- 
tion, ^B is the Bohr magneton, and j represents the jth 



First term (i^o) in the square bracket in Eqn.® is in- 
dependent of magnetic field, as because, under suitable 
variable transformation: /i — 2^bBx = y, we can elimi- 
nate B from i7o = 2^bB f{iJ — 2^BBx)dx^ and get 
the form rip = jl^o ^^v)dy- ^O'^' expanding the loga- 
rithm in f{y), and integrating over pz and j/, we recast 
f2o as 
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where At — y^Trfi^/mfc^jT is the thermal de Broglie 
wavelength, z — e}^I^BT |-j^g fugacity, and Licr(x) = 
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X + jl" + + ... is a poly log function of order 

a = 5/2. We already have mentioned, that, [i is the 
chemical potential of the system in absence of the mag- 
netic field. Thus, for a fixed average number (A^) of par- 
ticles (electrons), \j, is to be obtained using the implicit 
formula [2^ 




So far we have adopted the steps shown in the famous 
book of Landau and Lifshitz Q . Let us now proceed in 
our own way. 

B. Evaluation of integrals in the grand potential 



diamagnetism and Pauli paramagnetism 0, H, US] ■ This 
approximate formula, however, is not conclusive from 
Eqn.dH), as because, the approximation: hj ^ 1 although 
is applicable for lower values of j, but not for higher val- 
ues (> \) of j. Besides these, contributions of higher 
values of j are probabilistically negligible for negative 
values of the chemical potential, but not for positive val- 
ues. Thus, for 6 ^ 1, Equ.® is essentially true only for 
negative values of the chemical potential, but not conclu- 
sive for positive values of the chemical potential. 

On the other hand, for the case of very strong magnetic 
field (6 3> 1), we can use the expansion: 

k=\ 



The summation term in the square bracket in Eqn. (jS]) , 
on the other hand, can be evaluated expanding the log- 
arithm in /(/i — 2^,bBx) and integrating over pz and x 



(using the formula e cos(6a::)da 



to yield 
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where b = is a scaled magnetic field. The second 

KB-L ..... 

line of this equation is obtained inserting the constraint 
in Eqn.® and summing over k using the Mittag-Leffer 
expansion of coth(a;): ^ -l-2x ^^j^ -p^p^^. Since no ap- 
proximation has been made so far, Eqn.® represents the 
most general formula for magnetism of the system of our 
interest, and of course, is applicable for all possible values 
(strengths) of the magnetic field as well as for all tem- 
peratures. This equation must be the general formula for 
the Landau diamagnetism, Pauli paramagnetism and de 
Haas-van Alphen effect. An essentially similar equation, 
like this, was also obtained by Sondheimer and Wilson 
within a different (density matrix) formalism. But, this 
equation is difficult to work with for all temperatures and 
for all strengths of the magnetic field. While here-from 
they proceed for low temperature calculation [23| , we will 
proceed for all temperatures and fields. 

For the case of weak magnetic field (6 <C 1), Laurent 

series expansion formula of coth(j6) : ^ + ~ -^"^jg- + ■ • • + 

^^''(2k)\'jb ^ ■■■ frequently used, and Eqn.® is usually 

approximated within the lowest order in 6 as [6| 
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This approximate formula represents the magnetism of 
our system in weak field case, and unifies the Landau 



from the definition, and to the leading order in b, we 
approximate Eqn.® as 
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This equation leads to the saturation of magnetization of 
our system. 

Eqns.dTU]) and do not, of course, describe the mag- 
netism for all temperatures and fields, and are not true 
in general except in the two extreme cases. However, the 
expansion in Eqn. (1111) is true for all values of b. Using 
this expansion, and summing over j, we recast Eqn.® 
as 
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where v = fi/ksT. Ean. ((T5)) is also the most general 
equation like Eqn.® to describe the magnetism of our 
system. It is a novel form, and is easier to work with in 
comparison to the form in Eqn.®. In the next, we will 
evaluate the magnetization of our system for all values of 
magnetic field as well as for all values temperature using 
the form in Egn.dT^. 



C. Asymptotic expansion of the grand potential 

Sommerfeld asymptotic expansion of Fermi integrals 
(polylog functions) is of common interest for the low 
temperature properties of Fermi systems But, for 
the system of our interest, the Sommerfeld asymptotic 
expansion does not work well in particular for the oscil- 
latory part of the grand potential in Eqn. ([T^ . It works 
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only for v ^ oo resulting 
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where C(s, a) = X]^o ^/('^ + ''^Y ^ Hurwitz Zeta func- 
tion. It is easy to check, that, this asymptotic result dif- 
fers significantly from the exact result even for i/ 3> 1 . In 
the following, we will present a novel asymptotic expan- 
sion technique for the grand potential adoptingLandau's 
calculation for the de Haas- van Alphen effect [1, . 

It is nice to look at the b dependent part at the right 
hand side of Ean. ([T5|) . It appears as the same form as 
in the left hand side of Eqn.(IS|) which, of course, is a 
different form of Poisson summation formula. Using this 
summation formula in Egn. p^ and evaluating the inte- 
grations once by parts, we get 



^LiTpi) 
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Replacing the polylog function in Egn. dTSI) by its integral 
from, we get 



Im 



SNksTb 



T72)^2^io Jo 



Li3/2(-^)r( 



1 -(- Qy-i'^-'ixb) 



dydx 



(16) 



Now, replacing the integration variable a; by z = 2bx 
y — V and integrating over z, we recast Ean. (|16p as 
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where F(x) = e~^^ e^ dy is a Dawson integral. Our 
result in Ean. (IT9|) is not only a generalization of Pauli- 
Landau result, but also is an exact formula for T — )• 0. 
First two terms of this equation lead to Pauli-Landau 
result as well as Sondheimer- Wilson formula [1, [23j 
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The second term in Ean. ([T5)) as well as in Eqn. ipO)) rep- 
resents de Haas-van Alphen oscillations around the con- 
tribution of the first term which represents the com- 
bination of Pauli paramagnetism and Landau diamag- 
netism. It is clear from Egn. ipH)) . that the oscillatory 
term plays significant role only in the strong field regime 
(ksT < (ibB ^ fi). It dies out exponentially in the 
weak field regime (^^ ^1). This oscillatory term ac- 
tually was obtained by Landau in 1939 [1,[2^. This term 
forms the basis for the de Haas- van Alphen effect. In the 
low temperature regime {kgT <^ ji), Eqn. ipOl) . of course, 
works well for weak and strong field cases, but not for 
very strong field case {ksT <^ ^ < ^bB). On the other 
hand, our approximate result in Eqn. (jl9p works well for 
all strengths of magnetic field in low temperature, and is 
capable of predicting saturation of magnetization for very 
strong field case. Egn.dT^]) represent a novel asymptotic 
series expansion for the grand potential of the 3-D gas 
of electrons exposed in a constant magnetic field. Low 
temperature properties of this system can now be eas- 
ily obtained from this expression. In the following, we 
will compare Pauli-Landau result with ours evaluating 
magnetization of the ideal gas of electrons. 



III. EVALUATION OF MAGNETIZATION AND 
SUSCEPTIBILITY 



where the asymptotic series expansion of the incomplete 
beta function part, for large v, is given by 
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Using this expansion formula and evaluating the integra- 
tions for < y < & 1, we recast Eqn. (fTTl) as 
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Magnetization (Mk) of our system is to be obtained us- 
ing the definition M = —y^g-- Similarly, Pauli-Landau 
result on the magnetization (Mp-^fc) is to be obtained 
from Eqn.dini) as @ 
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Right side of Eqn. ipTj) has already been plotted in FIG. 
1 using an almost exact temperature dependent formula 
of the chemical potential [1^ . We will illustrate the tem- 
perature dependence of the chemical potential later in 
this section. Before that, let us see how the exact form 
of the magnetization looks like. 
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A. Temperature and field dependences of the 
magnetization 

Now, we obtain value (M = M{T,B)) of the magne- 
tization taking the partial derivative of 5fl in Eqn. (jl3p 

as 



M(T, B) 



Li|(-^) 
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where n = N/V is the fixed average number density of 
particles (electrons). It is clear from Egn. ip^ . that, as 
& — >■ cxD, only the first term in the square bracket con- 
tributes to yield the saturation value of the magnetiza- 
tion as Ms = n^B- We recast Eqn. ((22|) with this satura- 
tion value as 



M{T,B) = Ms 
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Li3(-e-) 
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This is the most general formula for the magnetization 
of our system. For h <^ 1, using Euler-Maclaurin sum- 
mation formula in Eqn. (j23p . we approximate M as 
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Although this approximate formula, to the first order 
in 6, gives rise to the unified formula of magnetization 
(Mp_|_d) for the weak field Pauli paramagnetism and Lan- 
dau diamagnetism [1, 0, : 



Mp+d{T,B) = 



2Ms hbBUi/2{-z) 
3 fcsT Li3/2(-2:)' 



(25) 



the third and higher order terms in b individually go to 
infinite for z > 1 with alternative signs. Thus, it is dif- 
ficult to judge whether Eqn. (P5| be the leading term for 
6 <C 1 until we compare Mp+d and M plotting in a figure. 
To plot the same with respect to the magnetic field for 
a fixed temperature, we need to know the temperature 
dependence of the fugacity (z) or the chemical potential 



B. Temperature dependence of tlie chemical 
potential 

Chemical potential of our system by no means can 
be obtained in an exact temperature dependent for- 
mula. But, an almost exact formula, for the same, using 
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FIG. 2: Thick and thin solid lines in (a) represent our re- 
sults (Eqn.((23l)) for T/Tf = 0.05 and 1.2 respectively. Thick, 
semithick, and thin solid lines in (b) represent Ean. (|23p for 
HbB /kBTp = 0.2, 0.3 and 1.03 respectively. Dotted lines 
represent respective weak field results (Eqn.((25])). 



Eqn.®, was previously obtained (within an approxima- 
tion technique) as [1^ 

/i(t) w ^l<{t)e{Q.^2■i -t) + n>{t)e{t - 0.723) (26) 

where t — T /Tp is our new scaled temperature, Tp — 
2mkB (37'"^»')^^'^ is the Fermi temperature of our system in 
absence of the magnetic field, 9{t) is a unit step function. 
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is the chemical potential (in units of ksTp) in the quan- 
tum regime. 
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is the chemical potential (in units of kBTp) in 
the semi classical regime, and Jz{t) = [I927ri3 - 

[36\/6 - 9\/2]7rii9/2 -f ^/l{l22?>S^TH^ + [1152\/2 - 

4608V6]^5/2ti5/2 _ [864 - m2^]T:H^Y'^]''\ 



Plotting of magnetization and comparison with 
Pauli-Landau result 



Introducing the above temperature dependence of the 
chemical potential in Eqn. (f26|) . in FIG. 2 we plot the right 
hand sides of Eqn.([23|). Thick and thin solid lines in FIG 
2 (a) represent our results fEgn. ip^ ) on magnetic field 
dependence of the magnetization of the 3-D ideal gas of 
electrons for T/Tp — 0.05 and 1.2 respectively. The thick 
and thin dashed lines in this figure represents weak field 
results (Eqn.din])) for T/Tp = 0.05 and 1.2 respectively. 
The oscillations of the magnetization around the respec- 
tive dashed lines represent de Haas- van Alphen effect. 
The oscillations, according to Eqn. (|23|) (and Eqn. (f20| as 
well) are expected to be dying out exponentially as the 
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Comparison for T/7V=0.1 

M(T,B)/M, 




Susceptibility (a): for 1/7^=0.1, 1.2 (b): for /Js^/knTf^Q.!, 1.03 
Tf Tf 

;r(T,B)— ;i'(T,B)— 
c c 




FIG. 4: Thick and thin sohd lines in (a) represent our results 
(Eqn.((29ll) for T/Tf = 0.1 and 1.2 respectively. Thick and 
thin solid lines in (b) represent EQn. (|29|) for fiBB /ksTp — 0.2 
and 1.03 respectively. The dotted line represents weak field 
limit (Ean. (|30[) V Here C = fionfi%/kB is the Curie constant. 



FIG. 3: Solid (our result), dotted (Pauli-Landau result) and 
dashed (weak field result) lines represent Ean. (|23p . Ean. pi|l . 
and Eqn.((25]) for T = O.ITf respectively. 



our system as 



chemical potential decreases with the increase of temper- 
ature. And, it is the case for the thin solid line. Thick, 
semithick, and thin solid lines in FIG. 2 (b) represent 
our results (Eqn. (|23l) ) on temperature dependence of the 
magnetization for ijlbB /ksTF = 0.2, 0.3 and 1.03 respec- 
tively. Dotted lines in this figure represent corresponding 
weak field results. 

For consistency, our results must have to match well 
with Pauli-Landau result, in particular, for the weak 
{^lBB <C ksTp) and strong field cases (fc^T < ^bB <C 
knTp) in the low temperature regime. To show how 
Pauli-Landau result differs from ours in the very strong 
field case {fJ-sB ^ kgTp), we compare field dependence 
of both the resuhs in FIG. 3 for T/Tp = 0.1. From 
this figure one can conclude, that, Pauli-Landau result is 
reasonably good even in room temperature for weak and 
strong magnetic field cases, but not for very strong mag- 
netic field case. Unlike Pauli-Landau result, ours have 
the saturation limit as well as the classical limit for the 
magnetization. 

While for T > Tp, the magnetization in FIG 2 ap- 
proaches the classical limit; for T < Tp, it either os- 
cillates around the weak field result or approaches the 
saturation limit. But, these oscillations are not clearly 
apparent from the figure, in particular, for the lower mag- 
netic field regime. To make it apparent, we can plot its 
derivative (susceptibility) . 



D. Susceptibility 



Susceptibility (x) of our system is to be obtained us- 
ing the definition x{T,B) = Mo ^^^qb^^ ■ Now, from 
Eqn. d^^ . we get the exact form of the susceptibility of 



X(T,B) = ^io 



k^ 



Li 



i3(-e'"^ ^ ^ ^ 



fc=i 



-kLU-e'^-''''^)} 
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This formula of susceptibility is valid for all strengths 
of magnetic field as well as for all temperatures. On 
the other hand, weak field limit for the susceptibility 

(xp+d(r,i?) 

as 



^^aWp^^TVB)^ is obtained from Eqn.^ 



V (T m - n Lii/2(-^) 

Xp+d(J , B) - fiQ — - 

SkB-l Ll3/2{~Z) 
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Introducing the temperature dependence of the chemi- 
cal potential in Eqn. ((26l) . we plot the right hand sides of 
Egns. iP^ . Thick and thin solid lines in FIG 4 (a) repre- 
sent our results (Ean. ((29)) ) on magnetic field dependence 
of the susceptibility for T/Tp = 0.1 and 1.2 respectively. 
Thick and thin dotted lines in this figure represents weak 
field limits (Eqn.dSOl)) for T/Tp = 0.1 and 1.2 respec- 
tively. The oscillations of the susceptibility around the 
respective dotted lines represent de Haas- van Alphen ef- 
fect. Thick and thin solid lines in FIG. 4 (b) represent 
our results (Eqn. (j29p ) on temperature dependence of the 
susceptibility for ^bB /kBTp — 0.2 and 1.03 respectively. 
Dotted line in this figure represents the weak field limit 
(Eqn. ([30| ). It is easy to check, that, for T — )■ 0, the sus- 
ceptibility does not only oscillates rapidly with l/B, but 
also diverges as 1 jB^^"^ . 

For the oscillatory nature of the magnetization, the 
susceptibility may have negative values which are evi- 
dent in FIG. 4. So, a natural question arises, whether 
the specific heat of our system may have negative values. 
To know the same we have to evaluate temperature de- 
pendence of the thermodynamic energy of our system for 
different values of the magnetic field. 
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IV. THERMODYNAMIC ENERGY 

Although low temperature specific heat of the system 
of our interest was calculated even for the interacting 
case [13, nil , yet the same for the ideal case has not sur- 
prisingly been calculated for all temperatures and for all 
fields as well. In the following we will also not calculate 
the specific heat, but, will calculate the thermodynamic 
energy of the system for all temperature as well as for 
all strengths of the magnetic field. Specific heat can, of 
course, be easily apparent from the thermodynamic en- 
ergy. 

If we do not want to restrict average number of parti- 
cles in Ean. (|T3| . then, this equation can be recast along 
with Eqn.® as 



6n 



2VkBT 



2mbS 1 

-— -Ll3 — Z) 



k=l 



Li3(- 



-Lii(-z) 



(31) 



Using the thermodynamic relation: SE 



6n 



T^^l^^y — M^^|T,y (for a grand canonical ensemble), 
we get magnetic field dependent part of the thermody- 
namic energy as 

oo 



SE = 



2VkBT 



{5[Li3(-z) + 2^Li|(-e''-2'=^)] 



fc=i 



-3Li|(-z)}+4fe2^fcLi.(-e''-^'=^) 



k=l 



(32) 



Now, the field dependent part of the thermodynamic en- 
ergy per particle {5U = SE/N) is given by 



SU 



knT 



{6[Li3(-z) + 2^Li3(-e''-2'='')] 



k=l 



Li3(-z) 
-3Li 5 {-z)}+ Ab^ ^ kUi {-e"" 



-2kb\ 



fc=l 



(33) 



Weak field expansion of the above can be obtained from 
the Euler-Maclaurin summation formula; and, it results 



5U 



keT 
Li3(-z) 



--Lu(-z) + -Li_3(~z) 



-— Li_.(-z)-H0(5«) 



(34) 



It is interesting to note, that, field dependent part of 
the total energy is the same as the field dependent part of 
the grand potential only in the weak field regime. Field 
independent part of the energy per particle, on the other 
hand, is given by [1, i, [H 



3 Li5(— z) 
2 Li3(-z) 



(35) 





Energy (a): forT/T-f =0.05, 0.4, 1.2 (b): for /ijjB/jtg 7^=0.05, 0.4, 1.03 

2. h 
1.5 - 

1 



FIG. 5: Solid, dotted, and dashed lines in (a) represent our 
results (Eqn.([36l)) for T/Tf = 0.05, 0.4 and 1.2 respectively. 
Solid, dotted, and dashed lines in (b) represent Eqn.(|36p for 
liBB/ksTF = 0.05, 0.4 and 1.03 respectively. 



Thus, total thermodynamic energy per particle of our 
system is given by 



U = Uo + SU. 



(36) 



We plot right hand side of Eqn. ([36|) in FIG. 5 to see 
the magnetic field as well as temperature dependences 
of the thermodynamic energy of the system. Although 
the energy has oscillations with respect to the magnetic 
field, it never decreases with the increase of temperature. 
Thus, specific heat of our system, unlike the suscepti- 
bility, never gets a negative sign. It is also interesting 
to note, that, energy weakly decreases in the weak field 
regime, oscillates with frequency ~ 1/B in the strong 
field regime, and increases linearly in very strong field 
regime. The last part of the interesting phenomena is a 
result of confinement of orbital motion of the electrons. 
It will not happen for spin magnets. 



V. CONCLUSIONS 

Pauli-Landau theory for the magnetism of the 3-D 
ideal gas of electrons is very much well established within 
the range of its applicability. But, it neither explains the 
saturation of magnetization, nor it does not have a clas- 
sical limit. So, it is not a complete theory. Goal of this 
paper is to give a complete theory for the same. Our the- 
ory unifies Pauli paramagnetism. Landau diamagnetism, 
and de Haas-van Alphen effect for the 3-D ideal gas of 
electrons in a single framework. It is applicable for all 
range of temperatures as well as for all strengths of the 
magnetic field. 

Our theory is consistent with the unified effect of Lan- 
dau diamagnetism and Pauli paramagnetism for weak 
field cases, and with the unified effect of Pauli paramag- 
netism, Landau diamagnetism, and de Haas-van Alphen 
effect for strong field cases. Our result goes beyond Pauli- 
Landau result, and is able to describe the range the de 
Haas-van Alphen oscillation within the entire range of 
the magnetic field. Our result is also able to describe 
quantum to classical crossover as the temperature is in- 
creased. 
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While Pauli-Landau result consists of oscillatory terms 
over the weak field limit, our result consists of oscillatory 
terms over the saturation limit. We also have done a 
novel asymptotic analysis of our result, in particular, for 
the low temperature regime. This asymptotic analysis 
gives rise to a generalization of Pauli-Landau result, and 
it eases comparison of Pauli-Landau (P-L) result with 
ours. Low temperature properties of the system can now 
be easily obtained from our asymptotic analysis. Ap- 
plicability of Pauli-Landau result or Sondheimer- Wilson 
formula, for a given strength of the magnetic field, can 
now be tested comparing with our result. 

The last section: thermodynamic energy, has been 
added to explore the thermodynamics of the ideal gas of 
electrons in presence of a constant magnetic field. From 
the expression of the energy per particle, one can now 
easily obtain specific heat, entropy, equation of state, etc. 

For electrons in a metal. Landau's formula (for the de 
Haas- van Alphen effect) gets a generalization to Lifshitz- 
Kosevich formula for the inclusion of the effective mass 
of an electron within the lattice structure of the ions in 



the metal 0, . Lifshitz-Kosevich formula is still used 
for the analysis of experimental data p^j - frzi [l9j even 
for the real (interacting Bloch) electrons in a metal with 
further modification due to interparticle interactions [sO] . 
But, like Landau's formula, Lifshitz-Kosevich formula is 
also not applicable for all temperatures as well as for all 
strengths of the magnetic field. Generalizing our result, 
for the interacting Bloch electrons in a metal, one can 
complete Lifshitz-Kosevich formula, and can apply the 
generalized complete formula for all temperatures and 
fields. 
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